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ABSTRACT 

Echelle spectroscopy has been carried out towards a sample region of the halo of the 
giant H ii region 30 Doradus in the Large Magellanic Cloud. This new kinematical data 
is the amongst the most sensitive yet obtained for this nebula and reveals a wealth of 
faint, complex high speed features. These are interpreted in terms of localised shells 
due to individual stellar winds and supernova explosions, and collections of discrete 
knots of emission that still retain the velocity pattern of the giant shells from which 
they fragmented. The high speed velocity features may trace the base of the superwind 
that emanates from the 30 Doradus starburst, distributed around the super star cluster 
R136. 

Key words: galaxies:starburst - galaxies: Magellanic Clouds - H ii regions - ISM: 
kinematics and dynamics - ISM:supernova remnants - ISM:individual (30 Doradus) 



1 INTRODUCTION 



The 30 Doradus nebula in the Large Magellanic Cloud 
(LMC) is the closest example of a giant extragalactic H ll 
region and the largest in the local group of galaxies. It is 
regarded as undergoing intense eno ugh star fo rmation to be 
referred to as a 'mini-starburst' by iLeithereg |1998:) and as 
such is an important nearby laboratory of both massive star 
formation and starburst phenomena. The high ly dynamic 
nebulosity (e.g. lMeaburnlll98ll : lMeaburnlll987l) is powered 
by a super star cluster of ~ 100 massive stars. Remarkable 
HST imagery of the environment of the ce ntral cluster of 
massiv e stars has recently been presented bv lWalb orn et alJ 
i2002h . The combined winds, UV radiation and supernova 
explosions from so many massive stars at a similar evolu- 
tionary epoch enables the generation of the nested giant (20 
- 300 PC diameter) shells that comp r ise the giant H ll re- 
gion jMeaburnlll980l : iMeaburiil Il990l : lLeithereilll99^ . On 
the largest scales, surrounding 30 Dor are supergiant (600 - 
1400 pc diameter) interstellar shells such as LMC3. 

The term 'shell' will be used in this paper rather than 
the commonly used term 'bubble' since it is preferable to 
use a term that is dynamically neutral and constrained to 
no specific geometry (e.g. spherical). The term 'bubble', of- 
ten erroneously presupposes a roughly spherical, pressure- 
driven, energy-conserving shell. This is certainly not the 
case for the supergiant shells which are unlikely to be either 
spherical or energy-conserving. The division between 'giant' 
and 'supergiant' when applied to the LMC shells will be for 
the diameter r anges above and re cently confirmed by the H i 
observations of lKim et al.l 1^9^. Different, though related, 



mechanisms must be involved in the formation of LMC shells 
in these distinctly separate diameter ranges. The most im- 
portant difference is that supergiant shells have diameters 
in excess of the neutral gas scale-height of the LMC. 

The overlapping giant shells comprising the halo of 
30 Doradus have been shown to be expanding at around 
50 km s"^ (e.g. [Meaburn 1984; IChu fc Kennicutd Il994l) 
whereas a multitude of 15 pc diameter regions exhibit out- 
flows of > 200 km s~^. The latter were interpreted as 
young supernov a remnants in the perimeters of giant shells 
llMeaburnl '1988^. Fig. 1 is a cartoon that illustrates the hi- 
erarchy of scale sizes present in a giant Hll region like 30 
Doradus. 

The brightest, dominant velocity components of 30 Do- 
radus are complex but seem to be comprised of three dis- 
tinct velocity regimes corresponding to H I sheets along 
the sightline. These are at 250 k m s~\ 270 km s~ ] and 
300 km s-i jMcGee et al.l 119781 IChu fc Kennicutd Il994l : 



iKim et al 3'1999^ In this paper, the systemic velocity, Vs 
taken to be the average heliocentric velocity (Vhel) of these 
components, 270 km s~^, in agreement with previous ob- 
servations (see for example Peck et al]ll997l : lMeaburnlll99ll : 



iGarav et aLlll993l : IClavtonl 



19871) . 



In this work, the aim is to investigate the faint high- 
est speed phenomena in the halo of 30 Doradus in order 
to complete the kinematical characterisation of this impor- 
tant giant Hll region. New echelle observations of the line 
profiles of the highest speed phenomena in the halo of 30 
Doradus have been made with unprecedented sensitivity. 
These are described in Section 2. The region of study was 
in the vicinity of the Wolf-Rayet star R130 and wEis se- 
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Figure 1. Cartoon of the 30 Doradus nebula to illustrate the hi- 
erarchy of shell scale sizes. The ambient density increases towards 
R136. 

lected as a representative portion of the halo of 30 Doradus. 
This area is located in 'sheU 3', of one of the several gi- 
ant shells that comprise the 30 Doradus region (n omencla- 
ture from 'Meaburn 198Cf: ' Chu fc KennicutlJll994l . see also 
[WaiiE & Hclfand 1991) and was o ne of the regions investi- 
gated bv lChu fc Kennicutd il994l) . Section 3 is a discussion 
of the high-speed motions and morphologies in terms of the 
structures and dynamics of the giant shells that comprise 
30 Doradus. How the high speed features relate to the out- 
flow of hot gas from the 30 Doradus nebula is also discussed. 
Conclusions are drawn in Section 4. 



2 OBSERVATIONS AND RESULTS 

A reproduction of the ESO PR Photo 14a/02 is shown in 
Fig. 2. This is an optical image of the 30 Doradus nebula 
obtained with the Wide Field Imager camera on the 2.2-m 
MPG/ESO telescope at the La Silla Observatory. A white 
box marks the region which is investigated spectroscopically 
in the present work. 

Spatially resolved, long-slit echelle spectra of 
[O III] 5007A emission lines were obtai ned using the 
Manc hester Echelle Spectrometer (MES) jMeaburn et alJ 
Il984^ at the 3.9M Anglo-Australian Telescope (AAT) on 
the night of 1994 October 17. The data were obtained 
using three parallel east-west slits simultaneously, sepa- 
rated by 45" which were then stepped southwards in 2" 
increments. The slit length was 165" and the slit width was 
300 fim{= 2" and 20 km s~^). The broad slit compromises 
spatial and spectral resolution but permits detection of the 
faintest, highest speed components of the line profiles. The 
exposure time for each multi-slit position was 1800s. The 
locations of all the slit positions are shown in Fig. 3 and 
the boundary of this figure corresponds to the white box in 
Fig. 2. In the discussion that follows, a slit is identified by 
which block of six it appears in. From north to south the 
blocks are labelled A, B, and C while within a block, also 




Figure 2. 30 Doradus nebula. This figure is a cropped reproduc- 
tion of ESO PR Photos 14a/02. The white box marks the region 
from which line profiles were obtained. The size of the region 
displayed is approximately 200 pc across 
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Figure 3. MES slit positions marked against the background 
nebulosity. This region corresponds to the white box in the pre- 
vious figure 



from north to south, the slits are numbered from 1 to 6. 
The data were processed in the standard manner. 

Greyscale representations of the position-velocity (pv) 
arrays of line profiles obtained from all the slits are displayed 
in Figs0J|S| These data are also presented and discussed in 
Al-Mostafa (1999). The scale is logarithmic in these three 
figures to enable all features to be discerned. In Figs |7| and 
|H| deep representations of the data from slits C4 and C6 
are displayed in order to highlight the faintest, high velocity 
material. 

These data are amongst the most sensitive obtained 
for the halo of 30 Doradus and reveal a wealth of faint. 
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complex high speed features. IChu fc Kennicut3 il994l) car- 
ried out echelle spectroscopy across the halo of 30 Doradus 
but many of the faint, highest speed phenomena revealed 
here were not detected. In the following discussion, general 
trends amongst the wealth of complex kinematical features 
are highlighted. Discussion of individual features may be 
found in Al-Mostafa (1999). Adopting a distance to the LMC 
of 55 Kpc gives that 1" « 0.27 pc. 

2.1 Bright systemic features 

The pv arrays are dominated by very bright continuous ve- 
locity features. These are due to the ionized overlapping H I 
sheets close to the systemic velocity being disturbed by the 
slowly expanding (50 km s~^) giant shells aro und 30 Do- 
radus dMeaburniligiol : IChu fc Kennicuttlll994 . These sys- 
temic kinematical structures have been investigated in detail 
in earlier work on the halo of 30 Doradus. 

2.2 Discrete high speed knots 

At the smallest scale (a few arcseconds or approximately 
one parsec), the pv arrays contain numerous localised high 
speed (±200 km s~^ with respect to the Vsys) velocity knots 
which do not appear to vary widely in spatial scale. They 
are visible in all the pv arrays and represent the finest-scale 
high speed substructure detected here. 

2.3 Velocity loops and arcs 

At a slightly larger scale (around ten arcseconds or a few 
parsecs), individual loops and arcs are discerned. They are 
not continuous in the pv arrays but are coherent velocity 
features made up of the discrete velocity knots. The clearest 
example is that seen at slit position C6 (Fig.|S|and the deep 
representation. Fig. (HJ. 

2.4 Large scale coherent velocity features 

At the largest scales (tens of arcseconds or approximately 
ten parsecs), high speed knots are found to trace out co- 
herent velocity features that slowly vary between being red 
and blue shifted. The clearest example is perhaps that in slit 
position B3 where at offsets of approximately to 50" the 
feature is redshifted and between around 50 and 100 " it is 
blueshifted with respect to the bright continuous feature. 

3 DISCUSSION 

3.1 Giant shells in 30 Doradus 

The current explanation for the origin of the 30 Doradus 
nebula is that as the massive stars within the nebula evolve, 
their winds (especially during the Wolf-Rayet phase) and 
subsequent s upernova explosions generate swept-up shells of 
ionized gas jMeaburni ri988. 19911). The shells are observed 
to have a hierarchy of sizes and velocities as one moves fur- 
ther into the halo of 30 Doradus. In the dense centre, the 
shell sizes and velocities are ~ 1 pc and ~ 10 — 50 km s~^ 
respectively, while in the halo the sizes reach ~ 100 pc with 
velocities of up to 100 km s~^ (see Fig. 1). The shells are 



prone to instabilities and can break up and fragment, vent- 
ing the interior pressure. For example a dense shell that is 
accelerating (due to either a rapid drop in the external den- 
sity or to a new supernova explosion within the shell) may 
break up via the Rayleigh- Taylor instability. Alternatively, 
dynamical overs tabilites can also lead to a s hell breaking up 
into fragments jMac Low fc Normanlll993h . In both cases, 
the fragments produced will have sizes of the order of the 
shell thickness. In general, the halo of 30 Doradus, into which 
the shells are expanding, is inhomogeneous. This inhomo- 
geneity, and the disruptive effects of nearby supernovae and 
winds from stars not within the original shell, will mean that 
a shell will not remain coherent for long. 

It is important to note that the LMC is thought to be 
flattened and viewed close to face-on. T he scale height of the 
H I in the LMC disk was calculated bv lKim et alJ llT999^ to 
be ~ 180 pc so that it is likely that the structure 30 Doradus 
may also be somewhat flattened. The scale-height imposes 
a limit on the sizes of the shells that can be formed, ir- 
respective of how intense and coeval in time the massive 
star activity is. As the shells grow, they be come elongated 
in th e direction of the density scale height JKoo fc McKed 
Il992i) . leading to a break up of the shell in this direction 
and a 'blow-out' of the hot interior gas into the galactic 
halo. The remaining structur e is known as a galactic chim- 
ney iNorman fc Ikeuchi||l989 |) and thev have been observed 
in the Galaxv llNormand^a^ ^^ .11999) and in the starburst 
galaxy M82 iWills et al,. 199^ . 

The giant shells of 30 Doradus are likely to be the 
maximum sized spherical momentum conserving shell struc- 
tures, since the scale height of the LMC is comparable to 
their diameters. The supergiant shells far exceed this scale 
height and may be collections of fossil chimneys viewed face- 
on and also the result of pr opagating star formation (see 
e.g. iMcCrav fc Kafatoslll987ri that is constrained to proceed 
in the plane of the galaxy, resulting in a ring shape super- 
giant shell. The loss of driving pressure means they are ex- 
panding in a momentum conserving phase and surround a 
low de nsity cavity. Such cavit ies are clearly seen in the H I 
data of iKim et all and lStavelev-Smith et all J2002l) . 

The hot gas that has escaped from the interior of the gi- 
ant shells will enter the LMC halo. There is stron g evidence 
for such a halo in the LMC. iWakker et al.l[r998f) have used 
GHRS/HST observations to detect C iv absorption towards 
LMC stars that do not to reside within a shell. This means 
the hot gas implied by these observations is n ot local to the 
star and is likely to reside in the halo (see also lSavaee et all 




3.2 Origin of the velocity features 

It would seem unlikely that the high speed knots represent 
random gas clouds within 30 Dor since that would require 
an explanation for both their hypersonic velocities and the 
systematic way they are distributed about the pv arrays. 
In terms of the scenario discussed above (section 3.1), a 
straight-forward interpretation of the kinematical features 
is as follows. The largest scale features represent old giant 
shells that have broken up via Rayleigh- Taylor (RT) insta- 
bilities. An instability is generated as the shell is accelerated 
by its interior pressure through a decreasing ambient den- 
sity. Those portions of the shell that are expanding in the 
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Figure 4. Position- velocity arrays of line profiles from slit block A. The vortical scale is heliocentric radial velocity, Vhel ( km s ^) and 
the horizontal scale is in arcseconds. At the distance of the LMC, 1" = 0.27 pc. 



plane of the LMC are less prone to disruption. The frag- 
ments that used to be part of the shell have continued to 
coast at the pre-break up velocity and together these re- 
main as a coherent velocity feature. For the RT instability, 
the characteristic knot size at the break up of the shell will 
be of the order of the thickness of the shell. The sizes im- 
plied by this picture seem reasonable - the old shell will have 
a dimension of up to a hundred parsecs towards the outer 
regions of the nebulosity while the shell wall will have had 
a thickness of a few parsecs. These estimates are in accord 
with measurements of intact shells observed wi thin the ha lo 
of 30 Doradus and elsewhere in the LMC (e.g. IOevlll99^ . 

The smaller chains of high-speed knots could be due 
to more localised disruptions of the giant shells due to, for 
example, a neighbouring supernova explosion. This latter 
scenario was proposed bv iR.edman et all (|l99flh to explain 
the unique Honeycomb nebula, which lies in the halo of 30 
Doradus. They argued that its cellular structure is due to a 



shell that has begun to fragment by a RT instability being 
impacted by a blast wave from a nearby SN explosion. There 
have been approximately ~ 40 supernova explosio ns within 
the halo of 30 Doradus in the last 10* yr alone l|MeaburDl 
;1991) (in comparison, in the starburst gala xy M82 there 
have been ~ 50 SNe in the last ~ 200 yr; iMuxlow et alJ 
119941) . 

In the H I pv array data of lStavelev-Smith et al.l i2002l) . 

the LMC and the Galaxy are well separated in velocity. In 
their data, the Galaxy does not exhibit kinematical features 
with a Vhel ^ 100 km s~^ while the LMC does not ex- 
hibit kinematical features with a Vhel ^ 100 km s~^. In 
our data faint velocity features are seen from the Vhel of 
30 Doradus down to a Vhel ^ 100 km s~^. However, it is 
unlikely that these velocity features are associated with the 
Galaxy rather than the LMC for several reasons. Firstly, 
such features are not seen at slit positions offset from 30 
Doradus; secondly, many of the features can be traced back 
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Figure 5. Position- velocity arrays of line profiles from slit block B.The vortical scale is Vhel ( km s ^) and the horizontal scale is in 
arcseconds. At the distance of the LMC, 1" = 0.27 pc. 



to the 30 Doradus systemic velocity; thirdly, there is no 
known Galactic H ll region or ionizing source along the line 
of sight that could excite the [O lll] 5007A emitting gas and 
there is also no extensive background [O ill] 5007A emission 
in the galaxy (compare Fi g. Q and |H] here with figure 9 of 
IStavelev-Smith et aljr2002l) . 



3.3 Superwind from 30 Doradus 

Starhurst activity can give rise to a 'superwind' due to 
the intense radiation fields, winds and supernova explosions 
caused by the star formation. The extensive high velocity 
features revealed here may be marking the very base of a 
superwind localised around the 30 Doradus complex. The 
escape velocity of gas from the LMC in the neighbourhood 
of 30 Doradus is around 150 km s~^ so that the high speed 
ionized gas from disrupted shells and giant shells is escaping 



the gravitational pull of 30 Doradus and is being ejected per- 
pendicularly to the plane of the LMC, along the line of sight. 
The ionization boundary due to the R130 cluster (perpen- 
dicular to the plane of the LMC) will depend on the distri- 
bution of the gas but an upper limit of a few hundred parsecs 
can be estimated by calculating the Stromgren radius due to 
an ionizing flux of ~ 10^^ s~^ from the ~ 100 O stars and a 
mean gas density of ~ 1 cm~^. Assuming an ejection speed 
of 200 km s~^, it will take of the order 2 x 10® yr to reach 
a distance of ~ 500 pc from the point of origin and thus 
escape the 30 Doradus region. The gas will rapidly recom- 
bine once it has passed the ioniziation boundary and will 
then not be visible on the [O ill] 5007A pv arrays. The dy- 
namical timescale of the remaining giant shell walls is much 
longer since their progress in the direction of the plane of 
the LMC is slower than the material ejected perpendicu- 
lar to the plane. High velocity H I clouds will be formed 
by the escaping material as it recombines and these clouds 
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Figure 6. Position-velocity arrays of line profiles from slit block C. 
arcseconds. At the distance of the LMC, 1" = 0.27 pc. 

may be detectable in high resolution and high sensitivity H I 
kinematical studies. Of course, the kinematics of gas ejected 
from the LMC rapidly beco mes highly complex due to th e 
interaction with the Galaxy iWakker fc van Woerderll997l) . 



4 CONCLUSIONS 

In this work the kinematics of sample region of the 30 Do- 
radus nebula have been investigated using the MES. This 
intensive study has revealed high speed velocity features 
throughout this region. Although the kinematics are com- 
plex, general patterns are discerned at three different spatial 
scales. Small coherent velocity features are present through- 
out the region. These knots are often found to form loops 
and chains in the pv arrays and at the largest scales, can 
form velocity features which vary slowly between red and 
blue-shifted emission. It is suggested that all of these fea- 
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The vertical scale is Vhel ( km s ^) and the horizontal scale is in 



tures are explicable in terms of the current understanding 
of the 30 Doradus nebula. Shells and giant shells formed 
by the winds and supernovae of massive stars form and are 
then disrupted in the energetic turbulent environment of the 
halo of 30 Doradus. The fragments of the shells retain the 
velocity pattern of the original shell and are observed as the 
small high speed knots. If this explanation is correct, then 
high velocity knots are likely to be found across much of the 
face of 30 Doradus wherever the size of the giant shells have 
exceeded the scale-height of the LMC and led to 'blow-out'. 
The whole 30 Doradus nebula is flattened and viewed face 
on. The high speed velocity fragments are likely to form the 
base of an outflowing superwind that is escaping the galaxy. 
This is a microcosm of the processes that are taking place 
in starburst galaxies such as M82 in which there are many 
super star clusters like 30 Doradus and whose combined out- 
put lead to the spectacular optical filaments that mark the 
M82 superwind. 
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Figure 7. Deep presentation of position-velocity arrays of line profiles from slit block C4 to highlight fainter features. 
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Figure 8. Deep presentation of position-velocity arrays of line profiles from slit position C6 to highlight fainter features. 
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